Fermionic superfluidity and spontaneous superflows in optical lattices 
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We study superfluidity of strongly repulsive fermionic atoms in optical lattices. The atoms 
are paired up through a correlated tunneling mechanism, which induces superfluidity when re- 
pulsive nearest-neighbor interactions are included in the Hubbard model. This paired superfluid is 
a metastable state which persists for a long time as the pair-broken process is severely suppressed. 
The mean-field phase diagram and low energy excitations are investigated in a square lattice system. 
Intriguingly, spontaneous superflows may appear in the ground state of a triangular optical lattice 
system due to antiferromagnetic frustration. 

PACS numbers: 03.75.Ss, 03.75.Lm, 03.75.Hh 



I. INTRODUCTION 

Recently, the strongly correlated properties of ultra- 
cold atoms in optical lattices have attracted great inter- 
ests in physicists The high tunability of the interac- 
tion strength between atoms, as well as the easy manipu- 
lation of the optical lattices makes the system realistically 
viable to implement a quantum simulator [4|. The quan- 
tum atomic gases in optical lattices make it possible to 
build a model system so that we can explore the corre- 
lated properties in many-body physics such as supercon- 
ductivity, quantum magnetism, quantum criticality, etc, 
and examine the related theoretical models 15^ 



A recent experiment [8| showed that a couple of 
strongly repulsive atoms occupying the same site of an 
optical lattice can be stabilized by damping the single 
particle tunneling. The lifetime of the pair increases sig- 
nificantly with the on-site repulsion U of the Hubbard 
model, which is quite intriguing as intuitively an attrac- 
tive force between particles is required to obtain a bound 
state. In the presence of a periodic spatial potential, the 
energy of a particle does not vary continuously but is re- 
stricted to particular ranges of values. A pair of strongly 
repelling particles can be stable because if it fell apart, 
the two isolated atoms would ensure kinetic energies that 
fall in a forbidden band|3|. Another experiment directly 
observed that for a pair of strongly repulsive atoms in the 
optical lattice, single-particle tunneling is severely sup- 
pressed by the requirement of energy conservation while 
atom-pair co-tunneling is permitted through the second- 
order quantum process'lO]. Although these experiments 
were performed on bosonic atoms, it is conceivable that 
they can be applied to fermionic atoms since no quantum 
statistics is involved. 

The formation of metastable atom pair with repul- 
sive interactions was first proposed by A.F. Andreev 
in his study of diffusion of impurities in quantum 
crystal[ll|, [l2|. The quantum liquid of repulsive parti- 
cle pairs in optical lattices has been discussed by several 
authors jl3Ml7l|. In a recent paper, Rosch et al studied 
metastable superfluidity of repulsive fermionic atoms in 
optical lattices jlSj]. Other authors attempted to explore 
the possibility of counterflow superfluidity and to con- 



trol spin exchange interaction of two-species ultracold 
atoms in a commensurate optical lattice fT9l - [2^ . In the 
fermionic superconductivity and superfluidity, a central 
concept is pairing. In Bardeen-Cooper-Schrieffer (BCS) 
theory, electrons pair up by an attractive interaction me- 
diated by phonons of the underlying crystals. The attrac- 
tion between ultracold fermionic atoms is provided by 
Feshbach resonances. Another pairing mechanism comes 
from correlated hopping, which occurs in fermionic tight- 
binding models. It has been suggested as a possible ex- 
planation for the high- Tc superconductivity |23l425l| . 

In this paper we explore the fermionic superfluid (SF) 
in a system consisting of repulsive atom-pairs in opti- 
cal lattices. We deal with a system partially filled by 
couples of spin-up and spin-down fermionic atoms. The 
lattice site is either occupied by two atoms or empty in 
the one-band Hubbard model. In the limit of the strong 
on-site repulsion U ^ t, the pair-breaking tunneling is 
suppressed as is revealed by the experiment [1). On the 
other hand, the atom-pairs may transport across the lat- 
tice through the second-order quantum transition As 
a result, the atoms always move in pairs. The effective 
Hamiltonian is obtained through the quantum perturba- 
tion theory, which is mapped to an anisotropic antifer- 
romagnetic (AF) model instead of the usual ferromag- 
netic model in the Bose-Hubbard model. The system ex- 
hibits superfluidity when the nearest-neighbor (NN) re- 
pulsive interactions are included. This paired superfluid 
is a metastable state because the pair-breaking process is 
severely suppressed. We investigate the mean-fleld (MF) 
phase diagram and low energy excitations for a square 
lattice system. It exhibits a gapless mode in the SF state 
and a gapful mode in the solid state. 

The AF feature of the effective Hamiltonian may lead 
to an interesting phenomenon when the underlying opti- 
cal lattice is triangular. As well-known, the ground state 
has a long-range 120° Neel order. The variation of the az- 
imuthal angles between the NN spins corresponds to the 
phase modulation of the superfluid state, which leads to 
a spontaneous superflow in the ground state due to AF 
frustration. As a result, the system exhibits a pattern of 
convection consisting of vortex-antivortex pairs. 

The paper is organized as follows: In Sec. II we obtain 
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an effective Hamiltonian for strong on-site interactions 
by the second order quantum perturbation approxima- 
tion. In Sec. Ill we demonstrate the superfliudity of the 
paired fermionic atoms by mapping the system to the 
pseudospin-1/2 antiferromagnetic model. The phase dia- 
gram and low energy excitations are calculated in the MF 
approximation. In Sec. IV we explore the possible phe- 
nomenon of spontaneous superflows of the paired atoms 
in a triangular lattice. A summary is included in Sec.V. 

II. EFFECTIVE HAMILTONIAN 

We focus on the partial pair-filling system with v < 1. 
For sufficient low temperature and strongly repulsive on- 
site interaction, the atoms will be confined to the lowest 
band which is described by the Hubbard model [5| , 

{ij)a i (ij) 

(1) 

where fia {flu) is the annihilation (creation) spin-cr 
fermionic operator, = n^^ + n^j^ is the number op- 
erator with riiu = fju fieri s-nd t^- is the tunneling matrix 
element. [/ > is the on-site interacting energy and 
1/ > is the NN interaction. In this work, we confine 
our discussions to the case oi U ^ tu,V. 

Since the pair-breaking processes are suppressed, the 
single particle hoping is eliminated in the second-order 
quantum perturbation theory. The on-site Hubbard term 



is considered as the the unperturbed Hamiltonain Hg. 
The hopping term is treated as the perturbation Hi, 
which should be calculated to the second order of ta/U 
to avoid pair-breaking. The NN interaction term com- 
mutes with the Hubbard term and will be included in 
the effective Hamiltonian later. Using a generalization of 
the Schriffer-Wolf transformation!^, 

H^Ho + + [S, Hi]] + ■■■, (2) 

where [S, Hq] = —Hi and S"^ = —S. Suppose |a) are the 
degenerate paired states of the unperturbed Hq with en- 
ergy Eq and 1/3) are the pair-breaking intermediate states 
of Hq with Ho\P) = El Then Ei=Eq- U. It should 
be emphasized that in the initial states |a) all atoms are 
paired up in the lattice sites while in the intermediate 
states 1/3) only one pair of atoms is breaking. We have 
{a\Hi\a') = (/3|i?i|/3') = and {a\S\l3) = {a\Hi\p)/U. 
Disregarding the intermediate states with more breaking 
atom-pairs which will involve higher order of ta/U , the 
the second order quantum perturbation Hamiltonian is 
then 

{a\H(^^\a')^^Y.^a\Hi\P){miW)- (3) 

An alternative method of the degenerate quantum per- 
turbation theory can be found in Refs.[l9l[27l|. 
The effective Hamiltonian is then. 



H, 



u 



T.^fyynf3i + f^^f^ifyh) + 



(ij) 



u 



U 



(4) 



where z is the number of the NN sites. In the second- 
order quantum perturbation, the intermediate virtual 
state 1/3) that breaks the atom-pair has a lower energy 
{—U) than the initial and final states |a), which induce 
an attractive NN interaction. In order to prevent the 
atom-pairs from congregation, a moderate repulsive NN 
interaction V is introduced in the original Hamiltonian 
([T]) to overcome the induced attractive interaction. The 
first term describes the pair-hopping, implying that the 
spin-up and spin-down atoms transport together across 
the lattice. The composite object behaves like a hardcore 
bosonic molecule because the fermion pairs always hop 



together to their nearest neighboring site and for each 
site only one pair of atoms is allowed. It should be noted 
that the pairing in our work is of s-wave type, where 
only the lowest single-particle band is considered in the 
Hubbard model. 

To study the MF properties of the system, it is con- 
venient to map the effective Hamiltonian Q to the 
spin representation [l^ by defining as Six — 
ifyit + /^t/a)/2, S,y = - f^tf^i)/2^, and 5,. = 

[Six,Siy]/i = (nj - l)/2. 



H 



8Mh 
u 



u 



E ^i^^i^- 

(ij) 



(5) 
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The NN interaction V also acts as an external magnetic 
field exerting on the pseudo-spins. 

In contrast to the ferromagnetic model [1^ [13] in the 
usual hardcore Bose-Hubbard model, the effective Hamil- 
tonian ([5]) represents an anisotropic AF model, where 
there are several competitive phases dependent on the 
pair- filling v as well as V and to-. Rosch et al ob- 
tained a ferromagnetic model by making a particle-hole 
transformation for the down spins [l^. At the MF level, 
we minimize the energy at fixed z-polarization or pair- 
filling. Suppose the classical spins are in the X-Z 
plane with an angle 9i to the z-axis, then a bipartite 
structure with sublattices A and B is employed to de- 
scribe the possible periodicity in the ground state [stI- 
[sPl. The candidate states include an easy-plane AF 
phase {9a — ^Ob) or paired SF with a non- vanishing 
order parameter {fi^fii) ^ 0, and a canted AF phase 
{cos 9 A ^ cos 9b), which is actually a checkboard solid 
with a non- vanishing in one sublattice while a 

vanishing (/if/ij,) in another sublattice. In addition, 
there is a phase separation (PS) regime caused by the at- 
tractive NN interactions. The easy-plane ferromagnetic 
phase {9a = 9b) is proved to have higher energy than the 
easy-plane AF phase and does not appear in this system. 



a usual checkboard solid instead of a supersolid. More 
accurate calculations such as quantum Monte Carlo sim- 
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FIG. 1: The V ~ u phase diagram for a square optical lattice 
system with = 0.1 and = l.ltj.. 



ulations have demonstrated that supersolid states indeed 
do not exist in a square lattice system [sills^. 



III. MEAN-FIELD RESULTS IN A SQUARE 
LATTICE 



We are now ready to explore the superfluidity of 
fermionic atoms in a square optical lattice (z = 4). Here- 
after we use the units oi U = 1. Figure 1 displays the 
V — V phase diagram for hopping integrals t^ = Q.l. We 
take, e.g., the ratio tt/^i = l-l- Other choice of t^/t^ 
does not alter the conclusion qualitatively. We compare 
the mean-field energies of each candidate phases to de- 
termine the ground state. Generally, the canted canted 
AF phase with cos 9a 7^ cos 9b 7^ 0, 1 is a supersolid. 
But in the square lattice we find 9b ~ or tt, imply- 
ing the supersolid order degenerates to an ordinary solid. 
We will reexamine this issue through the low energy ex- 
citations. The solid phase takes place in the regime of 
0.4 % V % 0.6. For the SF order, there is a 7r-phase 
difference between the two sub-lattices (canted AF or- 
der). The phase diagram is symmetrical with respect to 
the pair-filling v = 0.5 which results from particle-hole 
symmetry of the effective Hamiltonian ([4]) . 

We discuss the superfluidity in terms of the condensate 
density p,, = = \{h^f,C)?- The SF phase has a 

uniform condensate density ps = ^{l — v) for a given 
pair-filling, independent of the value of V . In Fig. 2, 
we plot the condensate density for the SF state (dashed 
curve) as well as the solid state (solid curves) versus the 
filling i/ for y = 0.1. At the MF level, the condensate 
density vanishes in one sublattice while does not vanish 
in another sublattice. This indicates that this phase is 
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FIG. 2: (Color online) Condensate density versus the pair- 
filling V for the superfluid state (dashed curve) and the solid 
state (solid curves). For the solid state, the two branches 
correspond to sublattices A and B. 

The superfluidity can also be investigated through the 
low-energy excitations, which provides an accurate probe 
for the nature of the quantum phase. We study the low- 
energy excitations by introducing pseudo-spin operators 
a\{ai) for sublattice A and b\{hi) for sublattice B, respec- 
tively. In this case, a\ = = f^^f-^, such as for sublat- 
tice A, is a composite bosonic operator. After making a 
rotation to align the local spins along the z-direction, we 
obtain the spin- wave type Hamiltonian in the momentum 
space as. 
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Hsw = '^{'^{2i cos 9 A cos 9 b + U sin 6*^ sin + ^k^-k + '^k^k + "k^p ^ ^7k("k^^k ^ '^k^-k ~ '^k^k ~ "k^k) 



(2t sin^A sin6'_B + U cos 9 a cos 9 B){a^a-^ + ^^^k) ;^(cos ^'Afljj.ak + cos ^sfoj^fej^)} + Hi 



linear : 



(6) 



r 



where = (cos k^+cos ky)/2, the renormahzed parame- 
ters f = 4^^^^/!/, U = A[V-{t^ + tl)/Ul = 2{zV-fi) 
and the summation for momentum k is restricted to a 
half of the Brillouin zone. Hunear includes the linear 
term in operators aj^ and b-^^. With Hunear = Q, 9a and 
9b are determined and the above MF result is recovered. 
This Hamiltonian ([6]) can be diagonalized in terms of the 
Bogoliubov transformation to obtain the low-energy ex- 
citation spectrum. For the SF state we have, 

w| - 4<(1 - 7j,)[i(l - + 2^.(1 - u){U + 2t)7j,]. (7) 



Figure 3 exhibits the excitation energies versus the mo- 
mentum along the a;-direction for ti = 0.1, = l.lij,, 
and V — 0.1. For the superfluid state in Fig. 3(a) with 
v = 0.2, the energy spectrum is linear at small momen- 
tum which reveals a gapless mode. At larger momenta, a 
energy dip plays the role of a roton as in the "^He super- 
fluid. The excitations of the solid state in Fig. 3(b) with 
1/ = 0.45 reveals a gapful mode. No gapless low energy 
excitation mode is found in this parameter regime. It jus- 
tifies that this phase is a usual checkboard solid rather 
than a supersolid. 
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FIG. 3: The excitation energies of (a) the superfluid state at 
1/ = 0.2 and (b) the solid state at u — 0.45 for the momentum 
along the x-direction. ti = 0.1, — l.lt^, and V — 0.1. 



IV. SPONTANEOUS SUPERFLOW IN A 
TRIANGULAR LATTICE 



We now study an intriguing phenomenon of sponta- 
neous superflow of the fermionic superfluidity for a tri- 
angular optical lattice system. Although there are some 
debates on the possible disordered ground state in a tri- 
angular AF model because of the geometric frustration 
as well as quantum fluctuations, the current consensus is 
that the ground state has a long-range Neel order (33l - [36| . 
We consider the classical spins Sj and employ the 120" 
XY-Neel order for three sublattices A, B, and C. Let 9a, 
9b and 9c be the corresponding polar angles which re- 
flect the spatial density variations |37H39| . the MF energy 
of the system is written as 



MF 



j-2 1 j-2 

-r 



) (cos 9a cos 9b + cos 9 b cos 9c + cos 9c cos 9a ) 



(8) 



2U 



(sin 6*^ sin 



sin 6*5 sin^c + sin^c sin6'A) + {V — —){cos9a + cos9b + cos^c), 

o 



where is the Lagrangian multiplier that controls the 
total pair-filling. In formula ([8]), the 27r/3 azimuthal an- 
gle difference between spins in the three sublattices has 
been incorporated. 

The MF energy should be minimized with respect to 
the angles 9 a, 9b, and 9c at a given polarization (Sz) = 
It has a form similar to that in Ref . ji^] except the 
first term may become negative. In that case, the system 



is phase separated. Generally, there is a supersolid phase 
with 9a ^ 9b ^ 9c- For V < + ti?/U, the ground- 
state is a uniform superfluid. We focus on the uniform 
superfluid phase {9a — 9b = 9c — cos~^{u — i)), which 
is implemented at a moderate NN interaction V ^ {tl^ + 
t'^)/U. We explore the implications of this 120° Neel state 
and its possible consequence in the paired superfluid. 
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orbital axis favors spontaneous superflows. Early the- 
oretical discussions of possible superflow in solid '^He in- 
volved quantum tunneling through ground-state vacan- 
cies, as well as Bose-Einstein condensation and quantum 
exchanges within the lattice J^]. In 2004, Kim 

and Chan reported the observation of the unusual su- 
perflow without resistance from frictional forces in crys- 
talline helium|44.. 45|]. This remarkable finding has now 
been confirmed [46l448l | . The latest experiments indicate 
that, rather than being an intrinsic property of a perfect 
quantum solid, superflows owe their existence to macro- 
scopic defects or extended disorder in the structure of 
solid helium. 

In a mismatched Josephson junction of ultracold 
fermionic atomic gases, M. L. Kulicfioj proposed an os- 
cillating superfiuid amplitude inside the weak link and as 
a result the so-called 7r-junction. If the junction is a part 
of the closed ring then spontaneous and dissipationless 
superfiuid current can fiow through the ring. 



SUMMARY AND DISCUSSIONS 



FIG. 4: (Color online) Spontaneous superflows of the paired 
superfiuid in the optical triangular lattice, (a) The azimuthal 
angles between spins correspond to the phase difference be- 
tween superfiuids on neighboring sites and a Josephson super- 
fiow spontaneously occurs along the edge, (b) The arrowed 
circles represent vortex- ant i vortex pairs which form the su- 
perflow convection. 



According to the spin mapping (/it/ij,) = {Si 



: sm W, e 



'i4>i 



the paired superfiuid has an order parame- 
ter phase as that of the azimuthal angle </>,; of the spin. 
Therefore, the 2tt/3 azimuthal angle difference of the 
spins implies a A(/)jj = 27r/3 phase difference between 
neighboring sites of the superfiuid. The Neel order of the 
antiferromagnetic model thus corresponds to a periodi- 
cal phase modulation in the superfiuid. In the theory of 
the Josephson tunneling, two weakly connected superfiu- 
ids or superconductors will induce a current as a result 
of their phase difference as J oc sin At/). Consequently, 
the ground state superfiuid spontaneously fiows along the 
link of the neighboring sites, as is shown in Fig. 4(a). In 
the triangular lattice system, the superfiuid fiows form 
a closed ring-like vortex. Figure 4(b) schematically dis- 
plays a regular convection pattern of superfiows which 
contains a sequence of the vortex-antivortex pairs. 

Similar cellular superfiows and periodic textures were 
suggested in the ^He — A superfiuid when a perpen- 
dicular magnetic field is applied to a sample slap[41|. 
The coupling between the superfiuid velocity and the 



We have studied the superfiuidity of strongly repul- 
sive fermionic atoms from a correlated pairing mecha- 
nism. The superfiuid is a metastable state with the op- 
tical lattice sites either doubly occupied or empty. The 
composite objects transport in the optical lattice through 
the second quantum processes via virtual pair-breaking 
states. It exhibits superfiuidity below a critical temper- 
ature. Phase diagrams and low-energy excitations in the 
square optical lattice system are investigated. Due to 
the AF frustration, the correlated pairs may result in 
an appealing spontaneous superfiow phenomenon in the 
triangular optical lattice system. 

Some authors explored the possibility of formation of 
the non-s-wave BEC through Feshbach resonance in a 
nonzero angular momentum channel on a lattice with 
double occupation j50- 52]. Varying the detuning of non- 
s-wave resonance can lead to various quantum phase 
transitions between the phases: S-wave BEC, non-S-wave 
BEC, conventional Mott insulator and orbital Mott in- 
sulator (with broken lattice symmetries). This becomes 
possible when the atoms are confined in the p-orbital 
Bloch band of an optical lattice rather than the usual s- 
orbital band. The new condensate simultaneously forms 
an order of transversely staggered orbital currents, rem- 
iniscent of orbital antiferromagnetism or d-density wave 
in correlated electronic systems but different in funda- 
mental ways. 

The NN interaction depends on the overlap of the Wan- 
nier functions between the NN sites. A moderate value 
of V ^ {t'^ + 1^) /U is sufficient to create the superfluid- 
ity in the system. An alternative way of generating NN 
interaction by the long-range dipolar interaction is also 
possible [H, 1131 • In order to detect the superfiuidity of 
the correlated pairs, photoassociation spectroscopy may 
be used 15 51. Interference of matter waves released from 
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the lattice has been used to probe the superfluidity of 
single atom condensation[56| . By tuning the interaction 
from repulsive into attractive, the fermionic atom pairs 
are converted into molecules. The sharp peaks will ap- 
pear in the interference pattern of the released bosonic 



molecules due to the presence of a SF fraction. 

This work is supported by the National Natural Sci- 
ence Foundation of China under grant No. 10874018 
and by the 973 Program Project under grant No. 
2009CB929101. 
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